YITP-02-48 



February 6, 2008 



Can a brane fluctuate freely? 

Akihiro Ishibashi and Takahiro Tanaka 
Yukawa Institute for Theoretical Physics, Kyoto University Kyoto 606-8502, Japan 

Abstract 

No, it cannot in the following sense if a self-gravitating vacuum brane is concerned. Once 
we write down the full set of linear perturbation equations of the system containing a self- 
gravitating brane, we will see that such a brane does not have its own dynamical degrees 
of freedom independent of those of gravitational waves which propagate in the surrounding 
spacetime. This statement seems to contradict with our intuition that a brane fluctuates 
freely on a given background spacetime in the lowest order approximation. Based on this 
intuition, we usually think that the dynamics of a brane can be approximately described by 
the equations derived from the Nambu-Goto action. In this paper we fill the gap residing 
between these two descriptions, showing that the dynamics of a self-gravitating brane is in 
fact similar to that described by a non-gravitating brane on a fixed background spacetime 
when the weak backreaction condition we propose in this paper is satisfied. 
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1 Introduction 



Branes, general relativistic objects extended in spacetime, have been discussed in various contexts. In String or 
M-theory Q, branes are thought to be a fundamental constituent element. In particular, D-branes on which the 
endpoints of open strings are restricted to lie have turned out to play an important role. In cosmology, topological 
defects such as domain walls and cosmic strings may have influenced the history of our universe significantly. 
Recently, so called braneworld scenario has attracted much attention, giving an intriguing perspective that our 
universe is embedded as a brane in a higher dimensional spacetime. In particular, in a model proposed by Randall 
and Sundrum ||, the self-gravity of a self-gravitating tension brane plays an important role to reproduce the 
conventional law of gravity approximately. 

The dynamics of branes in a curved spacetime has so far been investigated extensively. One handy approximation 
to describe the dynamics of a brane is just to ignore its self-gravity, and to treat it as a test membrane on a given 
background spacetime. As far as a vacuum brane is concerned, its equation of motion follows from the Nambu-Goto 
action. The brane is supposed to oscillate freely following this equation in the lowest order approximation in which 
the self-gravity is neglected. Then, turning on the effect of self-gravity, one evaluates the emission of gravitational 
waves, and takes into account the radiation damping effect as a correction. In the cosmological context this method 
has been frequently used for the estimate of the emission rate of radiations Ji|, |j| . 

One might anticipate that arbitrary accuracy will be obtained by the successive iteration of this scheme. One 
way to justify this procedure is to solve both the equation of motion for the brane and the bulk field equations 
simultaneously, and compare the result with that obtained by the successive approximation. Though this program 
is impossible to accomplish in general situations, there are actually several exactly solvable cases. Making use 
of the high symmetry of the background geometries, the full set of perturbation equations of a self-gravitating 
vacuum brane in 4-dimensional spacetime with infinitesimal width have been solved 0, |). The unperturbed 
configurations considered in these references are composed of maximally symmetric 4-dimensional spacetime and a 
vacuum de Sitter 2-brane. It has been shown that, once the gravitational backreaction is turned on, a vacuum brane 
loses its own dynamical degrees of freedom, and the perturbed motion is possible only while incidental gravitational 
waves come across the brane. Namely, the conclusion of those references was that a gravitating brane is not able 
to oscillate freely. On the other hand, when we consider a test non-gravitating brane, a small deformation of a 
brane in a maximally symmetric configuration is described as a scalar field on the unperturbed worldsheet, which 
obeys the Klein-Gordon equation with negative mass squared ||, 0. The degrees of freedom corresponding to this 
scalar field have not been identified in the former description in which the effect of self-gravity is fully taken into 
account. In this sense, the physical picture described by the dynamics of a self-gravitating brane is quite different 
from that by a non-gravitating one. 

In this paper, we shall revisit the perturbation dynamics of a brane coupled to the bulk gravity, aiming at 
filling the gap between these two ways of describing the brane dynamics. Just for simplicity, we shall restrict our 
attention to a vacuum brane of co-dimension one. 

In Sec. ^, we remind the description by means of a non-gravitating brane. In the first part of this paper (Sees. || 
- ^), after briefly summarizing the unperturbed background geometries in Sec. we provide the general formula for 
bulk and brane perturbations in Sees. ^ and ^, respectively. The formalism developed there is given in terms of the 
master scalar variable introduced in Refs. flit |l2f and its variant, which considerably simplifies the perturbation 
equations. In Sec. 5.1, we discuss the interaction between the motion of a gravitating brane and bulk metric 
perturbations, deriving the boundary conditions for the bulk perturbations along the brane from the perturbed 
junction conditions. We also derive the equation of motion for a gravitating brane, and compare it with the equation 
for a non-gravitating brane in Sec. 5.2. In the second part (Sees. |6| and ^), we address the issue that in which 
situation the test non-gravitating brane treatment can well approximate the dynamics of a self-gravitating brane. 



We propose a weak backreaction condition in Sec. 6.1. In the succeeding subsections, we examine perturbations 
of a de Sitter brane embedded in a de Sitter bulk spacetime (dS„+i -brane C dS„+2-bulk) to show that when the 
weak backreaction condition holds the dynamics of a self-gravitating brane reproduces the picture described by 
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a non-gravitating brane. In Sec. |7| we complete our explanation of the correspondence between the two different 
ways of the description of brane dynamics, by constructing the retarded Green's function for the perturbation of 
a self-gravitating brane, and by identifying a pole in the expression of the Green's function with the perturbations 
corresponding to brane fluctuations. We summarize our results in Sect. ||. In appendix B, we explicitly construct 
the global solutions for the perturbations in some exactly solvable cases dSa^-brane C dS^-bulk. 



2 Dynamics of a non-gravitating vacuum brane 

First, we recapitulate an approximate treatment given by considering a non-gravitating vacuum brane on a 
fixed curved background ]9|, |l^|. We are concerned with an n-brane E embedded in (n + 2)-dimensional space- 
time (M,cimn)- We use x M for the coordinates of the bulk spacetime, while x^ for the coordinates on the brane. 
The metric induced on the brane is given by q^ v = gMNd^x M d v x N . The brane E divides the ambient manifold M. 
into two parts, which we respectively call M.+ and M.-. Dynamics of a non-gravitating brane follows the equations 
of motion derived from Nambu-Goto type action 

S = -a [ d n+1 xJ- det(q^) + Ap f d n+2 x^/-det(g M N), (1) 

where a is the brane tension, and Ap is the difference of the values of the vacuum energy density between M.+ 
and M-. We assume that both a and Ap are constant throughout this paper. In the non-gravitating brane 
approximation, we neglect the effects of gravity caused by Ap as well as those by a. Starting with this action, and 
taking the variation with respect to a small deviation in the brane configuration x (x 11 ), we have the equation of 
motion for a non-gravitating brane as 

Ap 
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Here the extrinsic curvature is defined by K^ v :— — gLNd^x M d v x L Vjjn", with n N being a unit vector normal to 
E pointing toward M + , and Vm the covariant derivative with respect to the bulk metric ()mn- 
The perturbations of a brane configuration, x M — > x M + 5x M , can be decomposed as 

Sx M = Zf + Z ± n M , (3) 

where Z¥ is a vector and Z± is a scalar living on the unperturbed brane worldsheet. The deformations tangential 



K% = '-. (2) 



to the brane Z 1 ^ are nothing but worldsheet diffeomorphisms, and the only transverse deformation Z±n M is 
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physically relevant |)], |lC| ]. Namely, the perturbations of the brane configuration are described by a single scalar 
field Zi_. The deformation of the brane on a fixed bulk generates changes in the induced metric and also in the 
extrinsic curvature as 

= -2K^Z X , SK^ = V^„Z ± + (Rx^ v - K\K a ^j Z x , (4) 

where V M is the covariant derivative with respect to the unperturbed background induced metric q^u, and R± f _ l ± u := 
n n dpX d v x Rmlnk with R M lnk being the background bulk curvature tensor evaluated at the location of 
the brane. From the condition that SK^^ = 0, we obtain 

(□ (n+1) -m 2 )Z ± = 0, (5) 

where Ot n+ i) := V^V^, and the mass-squared is given by 

m 2 = -R XflXv q» v - K\K\ = , (6) 

n 

with being the scalar curvature of the unperturbed worldsheet. Thus the brane perturbation obeys Klein- 

Gordon type wave equation defined on the unperturbed worldsheet. For dS-branes, the behavior of the solutions 
of Eq. has been studied in detail in Ref . || . 
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3 Background geometry 



Next, we discuss the dynamics of a self-gravitating vacuum brane with co-dimension one. We assume that the bulk 
spacetime of the unperturbed background is maximally symmetric. Namely, it is given by [n + 2)-dimensional de 
Sitter (dS n +2), Minkowski (E„ + 2), or anti-de Sitter (adS„+2) space, depending on the value of the bulk cosmological 
constant A. As we will see, the unperturbed geometry induced on the brane also becomes maximally symmetric. 
Such solutions are well known, but we briefly recapitulate them in order to establish our notation. 

Anticipating future extension, we begin with more general framework merely assuming that the bulk metric 
admits the isometry of (not n + 2 but) n-dimensional maximally symmetric space K n with constant curvature K. 
The metric takes the form of 



dS (n+2) = ~9MNdZ M dZ N = g ab {y)dy a dy b + r 2 (y)da 2 (nK) , (7) 

where da 2 n K ^ = r )ijdx l dx : > is the metric of K n . We assume that the brane admits the same isometry as /C" so that 
the induced metric q^ takes the form of 

q^dx»dx v = a B 2 {-dr 2 + A 2 (t) da 2 {n K) } , (8) 

where r is a proper time normalized by a constant as, and aflA(r) is the restriction of the function r(y) on S, 
i.e., asA(r) = r(y)|s. We introduced as just for convenience, and we specify its value later. 

The Israel's junction condition implies that the difference of the values of the extrinsic curvature evaluated on 
both sides of S is related to the intrinsic energy-momentum tensor as 



[Knv] := K +fUI - K-pv = k 2 T^v Tg M „ ) , (9) 



1, 
n 

where k 2 is the (n + 2)-dimensional gravitational constant. This condition comes from the (p, ^-components of 
the Einstein equations. Hereafter, quantities with subscripts ± denote their value evaluated on the respective sides 
J\4±. We denote the difference of the values of a variable Q evaluated on both sides by [Q] := Q+ — Q- as used 
above, while the averaged value by 

Q-=\{Q+ + Q-)^ (10) 

On the assumption that the intrinsic energy-momentum tensor of a vacuum brane is given by T^ v = —crq^ Vl 
the junction condition (H) becomes 



Di r 



(11) 



where D± := n a D a and D a is the covariant derivative with respect to the orbit space metric g a b- In a similar 
manner we use an abbreviation D « to represent the covariant derivative projected in the direction tangential to 
the brane. Substituting Eq. (§) into the Hamiltonian and the momentum constrains, we also obtain 

T^K, V = ^y^M, (12) 
V„T% = 0, (13) 

where we assumed the (n + 2)-dimensional vacuum Einstein equations with cosmological constant, Gmn = Rmn — 
\RgMN = —\n(n + l)XgMN m the bulk. 

Now we focus on the case that Ai + and are both composed of locally maximally symmetric spacetime. In 
this case, the metric function r(y) satisfies 

- {D//r) 2 + (D ± r) 2 - K + \r 2 = , (14) 
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in respective sides, A4+ and M. Then, the equation of motion for the brane (11) reduces to 

K+{D„r) 2 _ n 4 (A + - A_) 2 + k s <j 4 + 2n 2 £ 4 a 2 (A + + A_) _ _rj_ 
r 2 An 2 k 4 a 2 a% ' 



(15) 



where rj takes a value —1,0, or 1. When r\ ^ 0, we choose the value of as so as to satisfy this condition. In the 
case of rj = 0, we simply set as = 1. 

Since r(y)|s = olbA{t) and D // = aB~ 1 (d T ) a D a on the brane, Eq. (|l5| ) reduces to K + (d T A) 2 — -qA 2 . We can 
immediately solve this equation to find that the brane has a locally maximally symmetric geometry. In all possible 
combinations of K and n, we can introduce Gaussian normal coordinates as 

ds 2 {n+2) = d X 2 + q tlv dx^dx v = d X 2 + a 2 { X ) {-dr 2 + A 2 (r)daf ntK) } , (16) 

for the bulk geometry, and thus 

r{y a ) = a(x)A{r) . (17) 

We specify the position of the brane at \ = xb, and therefore as — &(xb)- In these coordinates, the unit normal 
vector n a to the brane is naturally extended into the bulk by (d x ) a as the unit normal to x — const, hypersurface. 
Correspondingly, the extension of the unit tangent is given by a(x)~ 1 (d T ) a . From Eq. (jl4), we also have 

< i8 > 

or more explicitly (d x a) 2 = —\a 2 + ?/. 

The unperturbed background geometry is summarized as follows. 

• rj = —1; adS n+ i-brane C adS n+ 2-bulk: In this case, K must be —1 for the existence of a solution of 
Eq. @. The solution 

A(t)=cost, (for if = -1), (19) 

indicates that the brane is adS space. In order that Eq. ( Jig ) has a solution, A must be negative. For 
convenience, we introduce the bulk curvature radius I defined by |A| = i~ 2 . Then, the solution on each side 
M + or M- is given by 

aM^coshg) , (for A < 0) , (20) 
up to a shift of the origin of the x-coordinate. The embedding bulk is adS space. 

• rj = 0; E„ + i-brane C adS n+ 2-bulk or E n+ 2= There are two cases K = 0, and — 1. In both cases, the 
solutions 

A = constant, (for K = 0) , A(t) = r, (for K = -1) , (21) 
show that the brane geometry is given by Minkowski space. The warp factor is given by 

a( X ) = ^exp (|) , (for A < 0) , a( X ) = a B , (for A = 0) . (22) 

The bulk is therefore adS, or Minkowski space. 

• rj = 1; dSn+i-brane C dS n +2-, En+2-, or adS„+2-bulk: The scale factor on the brane is given by 

A(t) = sinhr , (for K = —1) , A{r) = expr , (for K = 0) , A(t) = coshr , (for K = 1) . (23) 

The warp factor is given by 

a(x) = ismh (|) , (for A < 0) , a( X ) = X , (for A = 0) , a( X ) = ^cos g) , (for A > 0) . (24) 

This means that de Sitter brane can be embedded into an ambient spacetime consisting of any combination 
of the two maximally symmetric spacetimes. 
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4 Bulk-perturbation 



In this section, we provide the bulk perturbation equations, following the gauge-invariant formalism for perturba- 
tions developed in Ref . . Perturbations of the present bulk geometry can be decomposed into the tensor-type, the 
vector-type, and the scalar-type components classified by the number of tensor indices tangent to the n-dimensional 
invariant space (/C n , da 2 n K i). The Einstein equations for each type of perturbations decouple from the others. We 
focus on the scalar-type perturbation since the displacement of a brane does not couple to the others as was shown 
in Ref. Q. We introduce the gauge-invariant variables for the scalar-type bulk perturbation. It turns out that 
these gauge-invariant variables can be described in terms of a single scalar master variable, and the perturbation 
equations reduce to a single equation for this master variable. 

The scalar harmonics S defined on QC n , da 2 n K s) is a function which satisfies 

(A + fc 2 )8 k = 0, (25) 

where A := DiD 1 , and £>i is the covariant derivative on [K n , da 2 n K ~ ) ). The scalar-type harmonic vector and tensor 
are constructed from §k as 

§ ks : = -rASk , §k« - T2 A-DjSk + -TySk • (26) 
The metric perturbations of scalar-type are expanded as 

Kb = fab^ > Ki = T fa S * * ' h V = 12 2r " ( H hij §k + H ) > ( 27 ) 

k k k 

where the coefficients f ab , fa, Hl, and Ht are functions on the 2-dimensional orbit space (M 2 , g a b)- In the following 
we abbreviate the indices k labeling the eigenvalues of the harmonics and the summation over them. 

Under the action of an infinitesimal coordinate transformation x M — > x M + £ M , the metric perturbations 
transform as 

h>MN — * hjtfN + Sh-MN = h,MN — Vm^N — V N^N ■ (28) 

The generator of infinitesimal gauge-transformation for the scalar-type perturbation is expanded as 

£ a =T a §, & = rLS;. (29) 

Then, the metric components, f a b, fau ^l, and Ht, transform as 

5 fab = -D a T b - D b T a , (30) 

6 f a = -rD a (-) + -T a , (31) 
\r J r 

k D a r 

SH L = -—L -T a , (32) 

nr r 

SH T = -L. (33) 
r 

From these transformation laws, we find gauge-invariant combinations of the variables as 

F:=H L + -H T + ^X a , F ab := f ab + D a X b + D b X a , (34) 
n r 

where a vector in AT 2 was introduced by 

T ( . r 
k 



•V„ : ~ [fa + l D aH T ) , (35) 



which transforms as SX a — T a . Note that Si and Sy vanish identically for k = mode, and also vanishes 
for k 2 = nK mode. The gauge invariant variables which contain expansion coefficients of the harmonics at these 
special values of k 2 do not exist from the beginning. Hence, the equations containing these variables are no longer 
valid. In what follows, we do not consider k 2 (k 2 — nK) = modes. 
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As shown in Ref. [|ll| , the bulk Einstein equations for the scalar- type perturbation reduce to a set of equations 
for the gauge-invariant variables F and F ab (See appendix A.). Among them, we can identify the constraint 
equations 

F\ + 2(n-2)F = 0, D b (r n ~ 2 F b a ) — 2D a (r n ~ 2 F) . (36) 

From these constraint equations, we can show that there exists a master variable in terms of which the gauge- 
invariant variables above can be expressed as fll|| 

2nr n - 3 F = |n + 2^£) c |fi, (37) 

< n r n, r 1 _ 
r n ~ 3 {F ab + 2(n-l)Fg ab } = < D a D b H —D b H —D a > £1 , (38) 



where □ := D a D a . From the evolution equations of the Einstein equations, the equation for the master variable Ct 
reduces to 

f D c r k 2 } ~ 

i □ - (n - 2) — D c - — \n = 0. (39) 
In the Gaussian normal coordinates ( |l6|) , the master equation above is expressed as 

A n - 2 d T A-^d T + ^+P 2 + ??^} T p (t) - , (40) 

a n ~ 3 d x a-^d x + ^(P 2 + W 2 ) } K P ( X ) - , (41) 

where we assumed a separable form of the solution as ft — T(t)1Z(x), and introduced p 2 as a separation constant. 
Here, we also defined 

n- 2 
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The bulk Einstein equations for scalar perturbations and the derivation of the master equation are briefly summa- 
rized in appendix A. 

5 Perturbation of a self-gravitating brane 

In this section, analyzing the junction condition, we discuss the interrelation between the bulk metric perturbations 
and the displacement of the brane. The junction condition provides the boundary condition for the bulk metric 
perturbations on the brane, which we need in finding the global solutions of the bulk perturbations. Simultaneously, 
we obtain an expression for the gauge-invariant brane displacement solely written in terms of the master variable 
from the junction condition. Hence, we find that there are no dynamical degrees of freedom for the brane fluc- 
tuations independent of the bulk metric perturbations. Combined with the Hamiltonian constraint, the junction 
condition also provides the equation for the displacement of a self-gravitating brane, which is to be compared with 
that for a non-gravitating brane obtained from the Nambu-Goto action. 

5.1 Coupling of bulk and brane-perturbations 

To treat perturbations of a self-gravitating brane in a gauge invariant manner, the bulk metric perturbations have 
to be taken into account as well as the perturbations of the brane configuration. The perturbation of the brane 
configuration, Zj_, transforms in the same way as a combination of the bulk metric perturbations X±_ under the 
infinitesimal coordinate-transformation (p9|); 8Z± = SXj_ = T±. Thus we immediately find a gauge- invariant 
variable corresponding to the displacement of the brane as 

Y ± := Z x - X x . (42) 
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Since Y/i := Z» — X// is of no physical relevance as mentioned before, we simply set Yn to zero in what follows. 

Taking the bulk perturbation into account, we find that the perturbations of the induced metric and the extrinsic 
curvature are given by 



Sq^u = -2K liV Z± + h^ v , 
SK^ = VJ7 v Zi_+(R ±tklJ ,-K a tl K a J)Z ± 



(43) 



(44) 



Since bK^ v is gauge invariant, we can express it solely in terms of the gauge-invariant variables F, F a0l and Y± 
as 0, 



5K T T 



{11 11 tp ^ 

- 2 D//F ±// + -n a D b F«» - -D±F a a + -K% F//// - D}Y ± + ^^f,) S , 



= kl-F 



rD, 



-D \F 



D„r 



-F\ 



(45) 
(46) 



ID, r 



2 r 



-F\ 



D„r 



-D 



II 



k 2 

nr 2 n(n + 1) 



k 



Y ± }SS l j + - s Y ± S l j . (47) 



Let us derive the boundary conditions for the master variable at the location of the brane. The perturbation 
of the induced metric can be made continuous across the brane 



[8q»u] = , 



(48) 



by appropriately choosing the coordinates tangential to the brane because the intrinsic geometry on both sides of 
the brane must be identical. From the junction condition (O), we have 



[6K» V ] = . 



(49) 



For further computation, it is convenient to use the gauge Z±_ = so that X± = —Y±_ on the brane. Then, from 
the metric continuity (Mq) we have 



[F] 



-Y\ 



and 



[F //// ]=-2[K\Y x ], 



(50) 



(51) 



where we have used K T T = a 1 (d T ) a D » n a . Since K T T = —r l D^r for the present background geometry, we find 



[F////+2F] =0. 

From Eq. (^) and the component proportional to S l j in Eq. (]47|), we have 

[^//J=o, 

[Y ± }=0. 

From the component proportional to SS i j in Eq. (|47|), we obtain 



where we have used the formula 



2 [D ± F] = 
D a D b r 



-F 



±± 



-A g a t , 



(52) 



(53) 
(54) 



(55) 



(56) 



and Eq. (|l8|). No new independent condition is obtained from Eq. (|45|) . Note that the results we obtained are 
gauge independent, even though they are derived by using a specific gauge. 
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Now let us express these conditions 
and (53). Firstly, we rewrite Eq. (B^) as 



and (]5q) in terms of the master variable f2, using Eqs. (p? 



,2 fc2 



= \Di-{n-2)^^D// + — \ fi +(n-l) 



Di 



(57) 



where we have used the formula (jl5|). Further, using the master equation ([Io|), we can write the condition above 



as 



(n-1) 



= — (P 2 + ^ 2 ) 



Secondly, from the condition (p3|), we have 



£>// (r 



= 0. 



Thirdly, we write down the condition (B5J) explicitly. From Eqs. (|37|) and (pq), we have 



2r n ~ 1 i^ 



r(D x r)D x - r(D/r)£>/ + — J- O , 



'{Ji-L + 2(n- = Di+2-^-D±U2 



Di r 



Applying D±_ to the first equation, we obtain 



2r r ' 



(n - 1) — —F + D±F 

r 



(D ± r) 2 D ± + r(D 2 ± r)D± + r(D ± r)Dl 

k 2 1 - 

-{D ± r)(D,r)D, - r{D/,r)D ± D„ + —D ± }fl . 



Using these relations, we can rewrite the condition (|55j) as 

k 2 — nK r 



Di n 



where we have used 



which follows from Eq. (fL8|) . 

Combining Eqs. (|9|) and (|62|), we find 



D±r\ D\r r/ 



D,0 



= 0. 



Thus, D±_fi is continuous across the brane. Then, from Eq. (pq), we have 

1 



- [p 2 + r)v 2 ) 
A* 



n 



-2(n-l)a B D±n, 



where we have introduced 



~2 

K CTOLb OtB 



(58) 



(59) 



(60) 



(61) 

(62) 
(63) 

(64) 

(65) 
(66) 



2n 2 

The Eqs. (|m|) and ( |65|) are the boundary conditions for the master variable imposed on the brane. 

We should note that the master variable Cl is not continuous at the location of the brane. We introduce a new 
master variable defined by 

W := a- (n -'^D ± n = T p (t)Q p ( X ) (67) 

where 

Q P ( X ) := a-( n -VD ± K p . (68) 
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From Eqs. ([Il]) and (|64|), we find that a quantity Q p defined by (^) is continuous, and so is W. In the bulk, Q p 
satisfies the equation 

-( n -^d x a n - l d x + — 2 (p 2 + W 2 ) \ Q p = . (69) 



From Eqs. @, @, and ©, we obtain 

(p 2 + W 2 ) [nl = -a n - 1 [D X W] . (70) 

Then, we can derive the junction conditions for W from Eqs. ( |64|) and (|65| ) as 

[W] - 0, (71) 

[D±W] = 2ijl ~ l \ w. (72) 
as 

Therefore W obeys 

£W := /a-C"" 1 )^^- 1 ^ - ^-H M <5( X _ _ _L (V^cM"^- 2 ^ + *^ | w = . (73) 

Although we exclusively use this variable W or i7 in this paper, it will be also worth noting that in the bulk 
$ := a~ 1 A~( n ~ 1 ' ) yV satisfies the equation for a massive scalar field, (□( n+2 ) — n\)& = with □(„+2) := V m Vm, 
and ^ :— a _ ^" _1 - ) 9 x (a™$) satisfies that for a massless scalar field, □(„+2) v I / = 0. Furthermore, the junction 
conditions for 'J are simply given by [iff] = [Dj_W] = 0. 

5.2 Equation of motion for Y± 

From the condition (|l2|), we can also derive an equation of motion for the brane fluctuations, which is analogous to 
that for a non-gravitating brane (|[). We will see, however, that the equation for the self-gravitating brane differs 
in that it has a source term written in terms of the bulk metric perturbations. 

Let us derive the equation of motion for the gauge invariant displacement Y x . From Eqs. (^) and ([!?]), the 
trace of 8K^ U is computed as 

8K\ = ({-D 2 n (Ef) D, - £ - m 2 ] Y x + j) § , (74) 



J := -y(n+z) D// (r"+ 2 F //x ) + U?f)F xx -D x F, (75) 



with 



where we have used the constraint equations (p6[). Then from the perturbation of the condition (]1 

5K\ = , (76) 

we obtain the equation of motion for Y x as 

(a (n+1) ~m 2 )Y x + J = 0, (77) 
with the mass-squared given by Eq. (|^). Since there is no jump in Y x as shown in Eq. (jHJ), we have replaced Y x 



with Y±. Similarly, we can show that the source term J is also continuous across the brane using Eqs. (53) and (55). 
Thus, we replace J with J. Alternative derivation of Eq. ((77|) is found in Ref. |l^| . The appearance of a source 
term J given by the bulk metric perturbations shows a clear difference from the equation for a non-gravitating 
brane (||). 



More strikingly, Y x can be expressed solely in terms of the master variable. Using Eq. (jllj), from Eqs. (50), 
(|5l|), and (|54|), we immediately have 

[Fan 2(n - 2)F] = -^^l~ K 2 aY x . (78) 
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Then, writing F// « and F in terms of f2, we obtain 

" * = - 4(n- a i)r»-3 ( D /l + ^T D II + ^ 

Hence, Y± is entirely expressed by the master variable f2. This equation indicates that the brane loses its own 
dynamical degrees of freedom, and a self-gravitating brane can oscillate only while bulk gravitational waves keep 
coming to it. Since this statement was first claimed in Ref. ||, there has been a worry about the validity of the 
approximate description considering a non-gravitating brane, in which we assume the existence of the dynamical 
degrees of freedom corresponding to wall fluctuations a priori. 



6 Description of brane dynamics in the weak backreaction limit 

In this section, we discuss the perturbative motion of a self-gravitating brane and compare it with that obtained 
by considering a non-gravitating brane, using the formulas developed in the previous sections. We first propose a 
condition so that the brane fluctuations can be described with neglecting the gravitational backreaction of the bulk 
perturbations. We call it the weak backreaction condition. Whether or not one can take the weak backreaction 
limit depends on the system that we are concerned with. Next, analyzing dS„+i-brane C dS„+2-bulk system for 
n > 3, we perturbatively obtain a global solution for the bulk metric perturbations which forms a discrete spectrum. 
We will see that the motion of a self-gravitating brane corresponding to this special mode reduces approximately 
to that of a non-gravitating brane. These observations will fill most of the gap between the two ways of describing 
the brane dynamics argued in the preceding section for n > 3. In the last subsection, we separately discuss the 
case with n = 2, i.e., the case with 4-dimensional bulk. 



6.1 Weak backreaction limit 

Toward a resolution of the paradox mentioned in the preceding section, here a comment on the gauge dependence 
is in order. We can see from Eq. ( [79|) that 

F abl F ~ 0(0) ~ 0(jiY±) . (80) 
Let us assume fi is small. Then, it also follows by definition that 

f a , f abl H L , H T ~ 0(X ± ). (81) 

When we choose the gauge such that Z± = 0, we have Y± = — X±, and therefore 

f a , f ab , H L , H T ~ 0(Y ± ) . (82) 

In this gauge the dynamical degrees of freedom of the brane fluctuations are completely transmuted into those of 
the bulk metric perturbations. On the other hand, if we choose the gauge such that X± = 0, wc find Y± = Z±, 
F ab = fab, fa = ~k- x rD a H T and F = H L + \H T . Therefore 

/a, fab, H L , H T ~ O (fl) - 0(jiY ± ) . (83) 

Under this gauge choice, the brane fluctuations are mainly described by Z±, and the perturbation of the bulk 
geometry stays small. This gauge choice provides a description relatively close to that obtained by considering a 
non-gravitating brane on a fixed background, if fi is sufficiently small. 

Based on the above observation, we propose that the condition that we can neglect the gravitational backreaction 
effects on the brane motion is specified by 

H < 1 • (84) 
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We note that Eq. ( |l5| ) can be expressed as 

4 = 1+ (!^W) 2 + ^, (85, 

where r\ = 1, 0, —1 represents the signature of the intrinsic curvature of the brane as before. 

For a E-brane, Eq. ( |85| ) becomes trivial if we choose rj — 0. But by considering the limit of E-brane from rj = 1 
or —1 case, one can see that /i 2 goes to infinity in the E-brane limit. Hence, the condition ( |S4| ) cannot hold for 
E-branc. For dS-, and adS-branes, the weak backreaction condition ( |S4|) holds when 

| [A] | > kV, or |A|»kV. (86) 

The junction condition also gives 



if", n 2 [A] 



(87) 



Thus, the former case in Eq. (86) is equivalent to the condition that the gap between the values of the extrinsic 
curvature on both sides of the brane is sufficiently small. 
One can also express Eq. ( |l5| ) as 

r, 1 , (l-n 2 [A]/(£ 4 a 2 )) 2 



« 2 A+ 2n 2 [A]/(fc 4 a 2 ) + 4n 2 A/(fc 4 a 2 ) ' 

Hence, apart from the exceptional case that |[A]| is extremely small compared with A , the weak backreaction 
condition (Q) implies that the curvature radius of the brane is much smaller than that of the bulk, i.e., £+ ^> a. 
In the exceptional case ( |A| ^> |[A]| ) only the second possibility in Eq. ( js"6| ) can be realized. In this case, we have 
a ~ £_ ~ £ + , and the brane becomes an almost totally geodesic hypersurface. When we impose Z2-symmetry 
across the brane, we have |[A]| = 0. Hence, under Z2-symmetry the weak backreaction condition is satisfied only 
through the second possibility in Eq. (|86|). 

6.2 Global solutions in the weak backreaction limit 

Here, focusing on dS„+i-brane C dS„+2-bulk system with n > 3, we solve the lowest eigenmode of the bulk 
perturbations assuming the weak backreaction condition (|84|) . We find that there is at least one discrete eigenmode 
in this system. In the succeeding subsection, we will estimate the amplitude of the source term J for this discrete 
mode, and we will show that the backreaction of the self-gravity on the brane motion is small. 
Just for convenience, we introduce a new coordinate £ in place of x by 

In terms of £, the warp factor is given by a = £cos(x/£) = ^/cosh£. Here, as we have used the same functional 
form of the warp factor on both sides M±, the \~ an d also ^-coordinates become discontinuous at the brane. 
Hence, in order to manifest this distinction, we associate + or — index with the coordinate x an d C if necessary. 
Using the new coordinate, we can rewrite Eq. ( |6^ ) as 

c cosh 2 c J V cosh cy 

The solution that is regular also at £ — > +oo is given by 

Q p = e 7r?,/2 r(l + ^)(coshCrPj(tanhC) . (91) 

Once we notice the relation P~ v '(tanh£) = {eT vlxl /T(l + ^))(2cosh(')~ I/ , it is easy to see that there is a solution 
regular at £ — ► — oo for p = iv in the absence of a brane. Since this solution is constant and hence nodeless, it 
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is the lowest eigenmode if it is normalizable. The standard normalization condition for discrete spectrum using 
the Hermitian measure is given by J d( {cosh £) ~ 2v Q p Q*t — 5 PP >. Hence, it is actually normalizable for n > 3. 
We separately discuss 4-dimensional bulk case (n = 2) in the succeeding subsection. We will anticipate that the 
eigenvalue of this mode will shift from p = iv once the gravity of the brane is turned on, but the corresponding 
mode itself will continue to exist. This shift of the eigenvalue can be calculated as follows. Setting p 2 = — v 2 + e, 
we seek such a shifted mode, and examine the behavior of Q p to the linear order in e. Let us introduce a function 
£(x)by 

(92) 



Q p = exp J d\' £(x'[ 



where £ is supposed not to be very large everywhere. Then, substituting this expression for Q p into Eq. fl69|), we 
have an equation for E to the lowest order in e as 



1 £) + l = 0. 



Integrating this, we find 



D±W {±) 



(±) 



= e£(±)(x±) = - £ 1 



K = 



±-kI±/2 



d X ' ± a"- 3 (x±), 



(93) 



(94) 



where £{±) and hence W(±) satisfies the normalizability condition at £ 
derivative of W at the location of the brane is evaluated as 



±oo. The jump in the logarithmic 



D±_W 

w 



t + - 2 



„n— 1 



d X ~a n ~ 3 (X-) + 
d X+ a n - 3 ' 

ttI + /2 

^iY, (for n>3) 



7T^ + /2 



d X+ a n - A { X+ ) 



\x + ) 



(95) 



where we have used the fact that £+ w or < in the weak backreaction limit to obtain the approximation 
in the second line. Then, from the junction condition ([72]), we have 



4 r(u- 



c = 



r (v) 

Hence, we find that the shifted eigenvalue is given by 



2v 



iv{l 



2^2' 



2 r(i 



/'-• 



as 



6.3 Correspondence between two different descriptions of brane dynamics 



(96) 



(97) 



Now we carefully examine the effective equation of motion for Y± given in Eq. (77) focusing on the specific mode 
obtained in the preceding subsection. Using the junction conditions for the master variable, we can express J given 
in Eq. (75) in terms of the master variable Cl. Inserting Eqs. ( |37| ) and ( |38| ) into Eq. (|75|), with the aid of Eqs. ( |6~I| ) 
and (^), we obtain after some calculation 

1 , 



J 



VP 



r]v 2 ) O 



n 



with 



O 



1 



. r 8 .r 8 

A -^-A 5 — 

OT OT 



4(n - l)a%A n - 3 

This expression can be also obtained by directly acting the operator 

k 2 



1*1 
n A 2 



K 
A 2 



Di 



D,r 



D 



II 



m 



(98) 
(99) 

(100) 
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on both sides of Eq. (|79|). By using Eqs. ( |7C| ) and (|72[), J can be also written by using W as 

J = — 2(n - l)a 2 gOW . (101) 



The equation (J^) involves the eigenvalue p 2 , which is determined only after solving the bulk metric pertur- 
bations. Hence, Eq. ( |77| ) cannot be regarded as the evolution equation for Y± in a usual sense. It will be rather 
appropriate to think that J is a measure of the deviation from the equation of motion for the non-gravitating brane 



case. For the specific mode whose eigenvalue is given by Eq. (97), we have the expression 



J={^ T -^(fy + { ^0[n], for n>3. (102) 

Since tl is of 0(fiY±) as discussed at Eq. j80|), we can conclude that 

J ~ 0((jY±) , for n>3, (103) 

for this discrete mode. Therefore when the weak backreaction condition \x <C 1 is satisfied, the source term 



J appearing in Eq. (77) can be regarded as a small correction to the equation of motion for a non-gravitating 
brane (||). 

In the above we have identified one discrete spectrum as the perturbation mode corresponding to the brane 
fluctuations. However, one may be still puzzled with the statement presented below Eq. (f79|); Even in the limit of 
the weak backreaction, "a self-gravitating brane can oscillate only while bulk gravitational waves keep coming to 
it." We will explain that there is a loophole in this statement, and will find that the picture given by considering a 
self-gravitating brane is consistent with that obtained by considering a non-gravitating brane. First of all, even if we 
start with a non-gravitating brane, we can easily imagine that the bulk gravitational waves are emitted as a result 
of oscillation of the brane at the next level of approximation. The amplitude of the emitted gravitational waves will 
be proportional to the brane tension and hence proportional to \i. As we have seen in Eq. (|8C|), the amplitude of 
gravitational waves is smaller by factor fj, compared with Yj_ in the description by a self-gravitating brane. Hence, 
there is no inconsistency as for the relative amplitudes between bulk gravitational waves and brane fluctuations. 
The remaining puzzle is the following. Although the brane described by a non-gravitating brane approximation 
can fluctuate without the incoming gravitational waves, it looks unable to fluctuate without incoming waves in the 
picture described by the self-gravitating brane. This puzzle is solved by showing that in fact there is a solution for 
the self-gravitating brane which satisfies no-incoming wave condition. 

To show the existence of a solution satisfying no-incoming wave condition, we solve the equation for the r- 
dependent part of the mode functions (MOT) with A = coshr, 



where L := L + v with L = 0, 1, 2, • • •. We consider the mode function that satisfies the positive frequency condition, 
(coshr) -1 ' T p (t) ~ e~ lpT , at r — > — oo. This mode function is given by 

(cosh t)-%(t) = e np Pi p (-t&nhT) 



e 



L 
-ipr 



F[ -L,L + l,-ip + l 



r(l-ip) V i + e- 2 - 

T(-ip)e- lpT ( ~ - 1 

V-^ [-L.L+ 1, ip + 1; =- 

T(~ip + L + l)T(-ip-L) V '1 + e 2 - 

T(ip) e lpT 



^F[-ip + L+l,-ip-L,-ip+l;- =- . (105) 

r(L+l)r(-L) V 1 + e 2 -' 

Then, the asymptotic behavior of the mode function is given by 

„-ipr 

(caBhr)-%(r) ^ ^—y 
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r^+oc T{-ip + L + l)T(-ip- L) r(i + i)r(-i)' 

Choosing the time dependent part in this way, the solution with the eigenvalue given in Eq. ( |97| ) asymptotically 
behaves as oc e~ lp ( T ~'^ at r — > — oo with \(\ ~ — r. Therefore the mode function becomes a function of r — |C| 
alone. This means that there are no-incoming waves. On the other hand, at r — > +00 there are two components 
which behave as oc e~ lp<T ~'^) and oc e lp ( T+ '^). Hence, outgoing radiation exists. 

Finally we point out that the eigenvalue given in Eq. ( |97| ) has imaginary part. This means that the amplitude of 
this solution decreases like exp[— 3(p — w)t\. This decrease in amplitude can be understood as radiative damping 
effect due to gravitational wave emission. 

6.4 4-dimensional bulk (n = 2) case 

For n = 2, there is no normalizable discrete spectrum of W corresponding to the brane fluctuations. But we can 
still construct a solution which satisfies the no-incoming wave condition. When we choose the r-dependent part 
of the mode function as given in Eq. ( |105| ), the no-incoming wave condition is given by (coshO^Qp = e ip KI at 
|£| — > 00. For n = 2, this means that the solution with the no- incoming wave condition becomes 

Q p = e ipC , forC>Cs, Q P = e~ ipC , for C < (b ■ (107) 



Then from the junction condition (72), we have 

p=-ifi. (108) 

For this eigenvalue, we find 

J - 0(Ai 2 l±) ■ (109) 

As in the cases with n> 3,Y± therefore approximately satisfies the same equation as in the case of a non-gravitating 
brane. The only difference from the cases with n > 3 is that the solution is excluded from a complete set of solutions 
for W since it is not normalizable. 



7 Green's function 

In the preceding section, we have seen that the solution corresponding to the brane fluctuations is not an eigenmode 
for the n — 2 case. Thus, one may suspect whether such a mode is physically relevant or not. In this section, we 
examine the late time behavior of the retarded Green's function for the master variable W. Recently, a similar 
technique was used in the context of braneworld cosmology ju], [l5| . We will see that the solution corresponding 
to the brane fluctuation discussed in the preceding section naturally dominates the late time behavior even if it is 
not contained in the eigenmodes. 

For simplicity, we concentrate on the dS-brane C dS-bulk case as before. The Green's function for each com- 
ponent in the harmonic expansion in the K = 1 closed chart satisfies 

C G(t, X ; r', X ') = ~S(r - r')S( X - X ') , (HO) 



where we substitute A = coshr and k 2 = L(L + n — 1) in the expression of C given in Eq. ([73]). 

To find the expression for the retarded Green's function, we first construct 5{r — r') by a superposition of the 
product of r-dependent part of the mode functions given in Eq. (|l05| ). From the expression in the second line of 



Eq. (105), we can see that this function does not have a pole with respect to p, although this fact is not manifest 
from the expression in the last line. For the purpose of constructing 5{r — r'), these mode functions are to be 
normalized with respect to the inner product 

r°° dr 

(^.^)=/_ M ^-^T Pl (r)T B (r). (Ill) 



15 



It is convenient to introduce T p (t) := T p (—t). In fact, the inner product between T p and T p i is calculated as 

(%,%,) = ^Kv-v') _ (n2) 

P ' V P T(-ip + L + l)T{-ip-L) 
Then, the delta function 6(t — r') is constructed as 

Here we have not taken into account the contribution from the discrete spectrum. For large negative t or large 
positive r', the integration contour can be closed in the upper half complex plane. In the upper half plane the 
integrand has poles at p — Li, (L — l)i, ■ ■ ■ ,i. The contribution from discrete spectrum should cancel that from 
these poles completely. Hence, the contribution from the discrete spectrum can be taken into account by shifting 
the integration contour so that it passes beyond p = Li. 

The Green's function can be constructed by assuming the following form 

G(r, X ; r', X ') = £ dp G P ( X , x')^?(-ip + L + l)T(-ip L)|^g> . (114 ) 

Then, G P (XiX') should satisfy 

a-^d x a^d x + \( P 2 + v 2 ) ) G P (X, X') = Six - x!) ■ (115) 



* 2 j 

We denote the homogeneous solution corresponding to the outgoing wave at £ — > ±00 by Qp (x), with £ given in 
Eq. (§|). Then we have 

G P (x,x') = ^ {Q P + \x)Q p -\x')Q{x- x') + Q P +) (x')Q p -Hx)0(x' - x)} , (ne) 

with 

w P ~ {d ( Q p -Hx)}Q P +) (x)-Q p - ) (x){d c Q P +) (x)}- (117) 

The zeros of W p occurs when Q p ~^ is identical to Q p + \ Namely, W p vanishes when the solution satisfying the 
outgoing boundary condition at £ — — 00 simultaneously satisfies the outgoing boundary condition at £ = 00. For 
n = 2, we found in Sec. OA that the solution with the outgoing boundary condition on both sides satisfies the 
junction condition when p = —i/j,. Hence there is a zero of W p at p = — i/i. For n > 3, when the imaginary part 
of p is positive, the outgoing boundary condition implies that the wave function takes the decaying component 
toward £ — ► ±00. Thus, the wave function W with the eigenvalue given in Eq. ( |97| ) satisfies the outgoing boundary 
condition. Therefore this eigenvalue corresponds to a zero point of W p for n > 3. 

Let us consider the situation that a certain event localized in spacetime excites the brane fluctuations and/or 
the gravitational waves near the brane. This situation will be heuristically described by adding a localized source 
term on the right hand side of Eq. ( |73] ) as C W = S. The solution is given by using the retarded Green's function 
as 

W = { d n+2 S:'G(x,S:')S(x') . (118) 



If we consider the behavior of the perturbations at a time sufficiently late after the event, we will be able to close 
the integration contour in the lower half of the complex p-plane. Then, the dominant contribution comes from 
the pole of Q v that has the largest imaginary part. The integrand on the right hand side of Eq. (116) has poles 
at p = Li, (L — l)i, • ■ ■ as well as at the zeros of W p . The poles at p = Li, (L — \)i, • • • do not give a dominant 
contribution at late time because of the exponential suppression of r-dependent part of the mode function for 
these special values of p. As is explained above, the zero of W p with the largest imaginary part is nothing but 
the eigenfunction corresponding to the brane fluctuations. Therefore the late time behavior is dominated by this 
eigenfunction. Although this eigenfunction itself has a diverging amplitude at |£| — > 00 for n = 2, the solution 
(118) is guaranteed to go to zero for large |£| because of the causality of the retarded Green's function. 
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8 Summary and discussion 



We have studied the dynamics of perturbations of a self-gravitating vacuum brane of co-dimension one, taking the 
(n + 2)-dimensional maximally symmetric configurations as the unperturbed background bulk geometry. In the 
first part of this paper (Sees. § - 1^), we have extended the general formulas for the bulk and brane perturbations in 
terms of the gauge-invariant master variable, including the boundary conditions for the bulk perturbations on the 
brane derived from the perturbed junction conditions. We have derived the equation of motion for a self-gravitating 
brane in more general cases than in literature, confirming that the derived equation is different from that for a 
non-gravitating brane derived from the Nambu-Goto action. The evolution equation for the brane fluctuations in 
the case of a self-gravitating brane does not close by itself. It has a source term J which is specified by the bulk 
metric perturbations. We have also shown that the perturbed brane motion Y± can be completely determined by 
the bulk metric perturbations. This fact shows that the brane does not have its own dynamical degrees of freedom, 
hence it cannot oscillate freely. The brane can oscillate only while incoming gravitational waves are present. These 
results are generalizations of the works Q |[ and they made it manifest that the way of describing dynamics 
of a self-gravitating vacuum brane looks crucially different from what one might naively expect from the analysis 
of a non-gravitating vacuum brane. 

In the second part (Sees. ^ and [?]), we have addressed a question, in which situation and in what manner these 
two ways of describing the brane fluctuations can be consistent, using the formulas developed in the first part. 
We first proposed the weak backreaction condition to clarify the situation in which the effect of the self-gravity of 
the brane to the bulk metric perturbations can be regarded as small, and hence the dynamics of a non-gravitating 
brane well approximates that of a self-gravitating one. The criteria for the weak backreaction condition we have 
proposed is that a non-dimensional parameter fj, proportional to the brane tension and the brane curvature radius, 
defined in Eq. ( |6rj ) , is sufficiently small. This condition requires that the brane radius is sufficiently small compared 
to the bulk curvature radius or that the brane is an almost totally geodesic surface with negligible difference in the 
vacuum energies on both sides of the brane. 

We have examined the perturbations of the dS-brane C dS-bulk system in the weak backreaction limit. We 
have shown that even when the amplitude of the incoming gravitational waves are of order 0(/i), the amplitude 
of the brane perturbation Y± can be O(l). Nevertheless a self-gravitating brane does not have its own dynamical 
degrees of freedom because the brane fluctuation can be entirely written by the bulk metric perturbations. As 
mentioned above, in literature we can find the argument that a self-gravitating brane can oscillate only while 
incoming gravitational waves arc present. 

We have pointed out a loophole in this argument. If we allow the amplitude of perturbations to increase 
unboundedly in the past direction, there is a solution which satisfies the no-incoming wave condition. The solution 
has a discrete spectrum for n > 3, while the solution in the case of n — 2 is not normalizable under the usual 
Hermitian measure. We have also shown that, when the weak backreaction condition holds, the gauge invariant 
brane fluctuations corresponding to this solution approximately satisfy the same equation as for a non-gravitating 
brane. The correction due to the effect of self-gravity becomes small of O(fiYx) for n > 3 and 0([i 2 Y±) for n = 2. 

Nevertheless, from the nature of the unbounded increase of its amplitude in the past direction, one might 
think that such a solution should be discarded as an unphysical one. For the case of n = 2, the solution is even 
unnormalizable. However, one can see that this increase of the amplitude in the past direction is a rather expected 
behavior. In the treatment of the self-gravitating brane, we solve metric perturbations and brane fluctuations 
simultaneously. Therefore the leading radiation reaction is automatically taken into account. Hence, the radiation 
damping reduces the amplitude of the perturbation in the future direction, which means that the amplitude should 
increase toward the past direction. 

To reinforce the justification of identifying this solution with the one corresponding to the brane fluctuations, 
we have considered a situation in which there is an external force which excites the brane fluctuations. Using the 
retarded Green's function for the perturbations of the self-gravitating brane system, we analyzed the perturbations 
induced by an event localized in spacetime. We find that the late time behavior after the event is dominated by 
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the solution with the no-incoming wave condition mentioned above. In this setup, this special solution becomes 
dominant only at a late epoch, and it is guaranteed that we do not see an unbounded increase of amplitude in the 
past direction because of the causality. 

In the present work, we have considered only maximally symmetric geometries as the background. It will 
be interesting to study less symmetric background case. Recently there has appeared an interesting work on 
the interaction between gravitational waves and cylindrically symmetric domain wall |ll| , which also provides an 
example in which dynamics of a self-gravitating brane can be well approximated by a non-gravitating brane. 

In the second part of this work, we examined only dS-brane C dS-bulk cases. An anti-de Sitter bulk spacetime 
plays an important role in the context of braneworld scenario. For the adS-brane c adS-bulk, and E-brane C adS- 
bulk cases, however, the spatial section of the bulk spacetime becomes open, and some subtle problems might 
appear in the boundary condition at the infinity. This issue will be discussed in a separate paper. Studies on the 
co-dimension two brane and the system composed of many branes will be also interesting. 
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Appendix 

A. The perturbed Einstein equations and the master equation 

The perturbed Einstein equations on a maximally symmetric spacetime can be reduced to the following four 
gauge-invariant equations on the 2-dimensional orbit space N 2 for modes with k 2 (k 2 — nK) ^ 0: 

-UF ah + D a D c F c b + D b D c F c a + {-D c F ab + D a F cb + D b F ca ) + (K ~ 2(n - 1)A^ F ab 

( D r Dhr \ 

-D a D b F c c - 2n D a D b + —D b + — —D a F 
\ r r J 

{ rr, D c r . , ,D c rD d r , , D c r , s k 2 - nK 
- I D c D d F + 2n—^-D d F cd + n(n - 1) F cd - 2nDF - 2n(n + 1 D C F + 2(n - 1) = F 

-DF C C - n^-D d F c c + ^F c c - (2n + 3)A J F C C | g ab = , (119) 

-^Db(r n - 2 F b a ) + rD a (Ja) + 2(n - l)D a F = , (120) 

1 Jn n LOf i\ Dar n , ( -iv ,D a rD b r D a D b r\ . 
-- < D a D b + 2(n - 1) D b + (n - l)(n - 2) h 2{n - 1) } F 

2 |_ r r z r J 

1 f_ , ,D d r n-lk 2 ,1 „ , , ( D c r (n-2)k 2 -nK] , s 

+- \ □ + (n - 1) D d g + A > F c + (n — 1) < □ + n D c - ± '- = >F = Q, (121) 

2 i r n r z J I r n ' J 

2(n-2)F + F a a = 0. (122) 



From Eqs. ( pCj ) and ( [L22[ ) we have 

D b {r n - 2 F b a )=2D a (r n - 2 F), F\ = -2(n - 2)F . (123) 
Then we can show that there is a master scalar ft = rQ, on M 2 by which F ab and F can be expressed as 

r"- 2 F= jL(D+2A)n, r n - 2 F ab =D a D b Sl- (?—±n + ?—^\)n, (124) 
2n \ n n I 
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which are equivalent to Eqs. ( p8| ) and (|37|). Substituting these into Eq. (119), we have 

(D a D b + Xg ab ) E(Sl) = , (125) 

where 

D-n D c - (n-2)\}Sl. (126) 

r r z J 

From this, we can finally obtain the master equation 

f)C 1.2 _ _ js } 

U-n D c = (n-2)X}n = 0, (127) 

r r z J 

which is equivalent to Eq. (|3^) . For more details of the derivation of the master equation, see Ref . (T^j . 
B. Perturbations of dS 3 5 -brane C dS 4 6 -bulk system 



Here as an exactly solvable system, we consider perturbations of 3, 5-dimensional de Sitter brane embedded in 4 
or 6-dimensional de Sitter bulk, i.e., n — 2 or 4, and A± = f ± 2 > 0. If we write down the radial equation by using 
the coordinate C defined in Eq. (|S9"|), we have 

dl+ P 2 + V { V -l)^{a- 1/ n p ) -0, (128) 

where v = (n — 2)/2 as used in the text. Without loss of generality we can set the values of (± on the brane to 
zero. Hence, for n = 2 or 4 the solution in the bulk is simply given by 

oT v Tl v = A±{p)e^ ± + B±{ P )e ±ipC± , (129) 

in A4±. We quote the junction conditions (Q) and ( |65| ) again as 

[d c n p ] = o, 

{ P 2 + v 2 )[R. p \ = -2(n-l)nd c K p . (130) 

dS3-brane C dS4-bulk case 

This system has been investigated in Refs. [Q, |). But for the convenience for the readers we quote the results here. 
We have from the boundary conditions ( |130| ), 

A = B, (131) 
^[A + B] = -A+ + B+. (132) 

Then, we obtain (see Eq. (4.49) in Ref. §) 

B± = -^A± + -^—A T . (133) 

P + l^L p + If! 

We see that for the pole mode p — — ifx, B± can be non zero even if A± = 0. Namely, there is a purely outgoing 
solution for p = —ifi. For this value of p, the source term is given by 

J = -U-rL- |5tanhr^- - L<yL + }" > + 4 tanh 2 T + u 2 lii, (134) 
^4a| \ dr 2cosh 2 r J 

where we have eliminated the second derivative of r, using the master equation. 
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dS5-brane C dS6-bulk case 



From the conditions (130), we have 



(p 2 -2)[A + B} + &ip^/l-\a 2 {A -B) + 3a 2 [\{A + B)] = , (135) 

and 

y/l - Xa 2 (A + 6)1 - 2ip(A - B) = . (136) 
From these, we obtain the relation between A±(p) and B±(p) as 

B± = UA± + VA T , (137) 

with 

/ill - 2p 2 + 3Jl - \+a 2 Jl - \_a 2 - Sip (jl - X+a 2 + Jl - \_a 2 ) } 

U := —± \ ^, (138) 

ip{p 2 + 1) - (4p 2 + l)/i - 6fi 2 ip - 3PV1 - \+a 2 ^Jl - A_a 2 

V := IP{P2 + 1) - (139) 

ip(p 2 + 1) - (4p 2 + l)/i - 6fi 2 ip - - A+aV 1 " A -« 2 

The zeros of the denominators are calculated to the linear order in fi as 

p=(ljA +a 2 A, - (l + ^A + a 2 M h, -{4-3A+a 2 }^, (140) 



where we have used 

y/l - \oA = [d x a] = -2fj, . (141) 
The eigenvalue with the largest imaginary part corresponds to that obtained in Eq. (^) . 
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